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ABSTRACT 

We present a complete spectral analysis of an XMM-Newton and Chandra campaign 
of the obscured AGN in NGC 4507, consisting of six observations spanning a period of 
six months, ranging from June 2010 to December 2010. We detect strong absorption 
variability on time scales between 1.5 and 4 months, suggesting that the obscuring 
material consists of gas clouds at parsec-scale distance. The lack of significant vari- 
ability on shorter time scales suggests that this event is not due to absorption by broad 
line region clouds, which was instead found in other studies of similar sources. This 
shows that a single, universal structure of the absorber (either BLR clouds, or the 
parsec-scale torus) is not enough to reproduce the observed complexity of the X-ray 
absorption features of this AGN. 
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1 INTRODUCTION 

The standard Unification Model for active galactic nu- 
clei (AGN) assumes the same internal st ructure for Seyfert 

2 and Seyfert 1 galaxies l|Antonuccilll993l ). with all the type 
1/2 observational differences ascribed to an axisymmetric 
absorber/reflector, located between the broad line region 
and the narrow line region, in order to obscure the former, 
but not the latter. An early developed, natural geometri- 
cal and physical s cenario is that of a homo geneous torus 
on a parsec scale (|Krolik fc BegelmanI flQSj '): however re- 
cent studies on X-ray absorbing column density changes per- 
formed with Chandra, XMM-Newton and Suzaku satellites 
ruled out a universal geometrical structure of the circum- 
nuclear absorber. Absorption variability is common (almost 
ubiquit ous) when we comp are observations months to years 
apart ijRisaliti et al.l l2002l ). and, most notably, has been 
found on time scal es of hours to days in several s ources, 
such as NGC1365 jRisaliti et"ai] |2005|. [20071. |2009|) NGC 
4388 (lElvis et all I2004I'). NGC415 1 (jPuccetti et all l2007l l 
and NGC 7582 (jBianchi et all |2009| ) . 

NGC 4507 is a nearby (z=0.0118) barred spiral galaxy 
and one of the X-ray brightest Compton-thin Seyfert 2s, de- 
spite the heavy obscuration (A'^e ~ 4 
first obs erved in the X-ray s by Einstein (Lx 
iKriss et all [l980l ) 
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^). It was 
2.8 X 10"^ 
and then in 1990 with Ginga 



(|Awaki et all [l99ll '). showing a strongly absorbed (A''h ~ 
5 X lO'^^cm"'^) power law continuum and a prominent iron 
Kq line. In 1994 ASCA also revealed a strong X-ray ex- 
cess and an intense emission lin e (identifi e d as Ne ix) at 
~ 0.9 keV JComastri et allll998l 'l. iRisalitil l|2002l 'l reported 
then three BeppoSAX observations of the source, confirm- 
ing the obscuring column densities observed in the previous 
two observations. In 2001 NGC 4507 has been observed for 
the first time with the most sensitive X-ra y satellites, XMM - 
Newton and Chandra/ACIS-S HETG |Matt et alj l2004t ) 
confirming once again the clear Compton-thin state of the 
spectrum, with a F = 1.8^0 2 power law absorbed by an 
JVh = 4.4t° j X lO^^cm"^. In 2007 a ~3 days observation 
with Suzaku has been performed and it revealed a much 
larger absorbing column density (A^h ~ 9 x lO^'^cm"^) with 
respect to the earlier ob servations, but no changes within 
the 3 days of monitoring (|Braito et al.ll20l3 . hereafter B12). 
Summarizing, NGC 4507 showed strong variations in time 
scales of years (from 1990 until 2007) but none in shorter 
time scales (3 days observation with Suzaku). This excluded 
a possible sub-parsec scale absorbing structure, but leaves 
the actual size and distance of the absorber largely uncon- 
strained. In this paper we report the study of five XMM- 
Newton observations spanning a period of 6 weeks, between 
June and August 2010, and a Chandra observation per- 
formed 4 months later, in December 2010. This observational 
campaign has been designed to fill the time gap between 
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days and years in the previous observations, so constraining 
the location of the absorber. 



2 OBSERVATIONS AND DATA REDUCTION 

The 5 XMM-Newton observations analysed in this paper 
were performed on 2010 June 24 (obsid 0653870201), July 

3 (obsid 0653870301), July 13 (obsid 0653870401), July 23 
(obsid 0653870501), August 3 (obsid 0653870601) with the 
EPIC CCD cameras, the PN llStriider et al■ll2001^ and the 
two MOS (|Turner et al.l [200 if) , operated in large window 
and medium filter mode. The extraction radii and the opti- 
mal time cuts for flaring particle ba ckground were computed 
with SAS 11 (jCabriel et al.ll2004h via an iterative process 
which leads to a maximization of t he Signal-to-Noise Rati o 
(SNR) , similarly to that described in lPiconcelli et all l|2004h . 
After this process, the net exposure times for the 5 different 
observations were 16 ks, 13 ks, 13 ks, 13 ks and 17 ks for the 
PN, respectively. The resulting optimal extraction radii are 
40 arcsec for the first three observations, 30 arcsec for the 
fourth one, 26 arcsec for the last one. The background spec- 
tra were extracted from source-free circular regions with a 
radius of about 50 arcsec for all the 5 observations. We also 
re-extracted the data from a previous XMM-Newton obser- 
vation (obsid 0006220201), with a net exposure of about 36 
ks, adopting an extraction radius of 40 arcsec for the source 
and 42 arcsec for the bac kground. The ana lysis of the last 
set of data is discussed in iMatt et all \2004 ). 

Chandra observed the source on December 2, 2010 for 44 ks, 
with the Advanced CCD Imaging Spectrometer (ACIS: 
ICarmire etallbOOSh . Data were reduced with the Chandra 
Interactive Analysis of Observations (CIAQ: rFruscione et al.l 
120061 ) 4.4 and the Chandra Calibration Data Base (CALDB) 
4.4.6 database, adopting standard procedures, using a 2 arc- 
sec and 10 arcsec extraction radii for the source and back- 
ground, respectively. 

Spectra were binned in order to over-sample the instrumen- 
tal resolution by at least a factor of 3 and to have no less than 
30 counts in each background-subtracted spectral channel. 
This allows the applicability of the statistics. 



3 DATA ANALYSIS 

The adopted cosmological parameters are Ho = 70 km s^^ 
Mpc~^, nA = 0.73 and Qm. = 0.27, i.e. the default ones in 
XSPEC 12.7.0 l|Arnaudll 19961 ). Errors correspond to the 90% 
confidence level for one interesting parameter (Ax^ ~ 2.7), 
if not otherwise stated. 



3.1 EPIC PN/MOS spectral analysis 

The soft 0.5-3.0 keV spectrum presents a strong 'soft ex- 
cess', which appears dominated by emission lines from an 
highly ion ized gas, as obser ved in most X-ray obs cured AGN 
l| Turner et al. 1997; Gua inazzi fc Bianchil |2007| ). Emission 
lines from H-like and He-like C, N, O, and Ne, as well as 
from the Fe L-shell, have been detect ed and reported in pre- 
vious obs ervations ([Matt et al.ll2004l ). 

Following iMatt et al.l l|2004h the baseline model we used to 
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Figure 1. XMM-Newton EPIC PN 0.5-10 keV best fit and resid- 
uals for tlie 5 observations of the campaign described in this work. 



fit the 0.5-10 keV spectra can be roughly expressed by the 
following general formula: 

F{E) = e"'^^^< [Phc + C + e-'^(^)'^«Bi5-^ + (1) 
+i?(r) + ^G.(£)] 

i 

wher e crjE) is the photoelectr ic cross-section (abundances 
as m lAnders fc Grevesselll989l ). A^^ is the Galactic absorb- 
ing column density along the line of sight to the source 
l|Dickev fc Lockmanlll99(il) : Phc is the emission from a pho- 
toionised gas repr oduced with self-con siste nt cloudy mod- 
els as described in Bianchi et al.1 (|2010l ) and lMarinucci et al.l 
(|201ll ) while C is the e mission from a col lisionally-ionised 
diffuse gas (apec model, [sinith et aLll200ll ): A^h is the neu- 
tral absorbing column density at the redshift of the source; 
B is the normalization of the primary powerlaw with slope 
F; -R(r) is the Compton scattering from the inner layer of 
the circumnuclear torus, mod elled in Xspec with pexrav 
l|Magdziarz fc Zdziarskilll995h :G, (i?) are Gaussian profiles, 
corresponding to required emission lines of high-Z elements 
such as the Fe Ka at 6.4 keV, Fe K/3 at 7.058 keV, Fe xxvi 
at 6.966 keV and the forbidden line of the Fe xxv Ka triplet 
(see table [J). 

A further Gaussian emission line has been used to reproduce 
the Compton Shoulder (CS) redwards of the Iron line core, 
as expected on the oretical grou nds, with energy fixed at 6.3 
keV and cr=40 eV (|Mattll2002l ). 

The reflected, Compton scattered emission has been mod- 
eled using the pexrav model with r=1.8 and normalization 
fixed to the values measured with the broadband Suzaku 
observation (B12). The previous model has been used to fit, 
in first place, the 5 separate EPIC PN/MOS observations. 
The 5 spectral fits from our campaign in Summer 2010 (la- 
beled as Obs. 1-5 in Table [T| have good x^/d.o.f. and do 
not present any strong evidence of variations in either the 
absorbing column density or F. The Obs. 2 data set presents 
strong residuals mainly between 1.5 and 3.0 keV, more ev- 
ident in the EPIC-MOS spectra. Since they are located in 
very narrow bins at ~ 2.8, 1.6 and 1.9 keV and they 
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cannot be ascribed to known spectral features we believe 
they are due to background/calibration issues. 

Two different photoionised phases and a collisional one 
are needed to model the 0.5-3.0 keV spectra of NGC 4507. 
Table [T] clearly shows that the soft X-ray emitting gas has 
not varied during the monitoring. This suggests that the 
interveining absorbing material, responsible for the column 
density variation between June and December 2010, might 
be much closer to the X-ray source than the circumnu- 
clear matte r responsible f o r the soft emission. As already 
discussed in lBianchi et al.l (|2006t ) and lBianchi fc GuainazzH 
|2003) this gas is likely coincident with the NLR. Further 
studies on the phoionisation mechanisms and extended emis- 
sion in NGC 4507 will be discussed in detail in the future 
(Wang et al., in preparation). 

We then analysed the 10 spectra (5 EPIC PN and 5 
MOSl-l-2, labeled as Set 1 in Table 1) simultaneously, using 
the model described above and linking all the parameters, 
except for the flux normalizations. The baseline model re- 
produces the 5 sets of data and some residuals are present 
around 1.8 keV. Indeed, the addition of a line at 1.77±0.03 
keV is required (Ax^ ~ 44, with a significance greater than 
99.99%, according to T^-teslQ); it can be identified as Si Ka, 
with a corresponding flux of 2.0 ± 0.4 x 10~® ph cm~^ s~^ 
. The primary powerlaw (F = I.SIq j) absorbed by a col- 
umn density of 9.0 ± 0.5 x 10'^'^ cm"'^. The photon index is 
in agreement with previous studies on this source, on the 
contrary a clear variation in the Nh can be noticed with 
respect to the old 2004 XMM-Newton observation, but it is 
fully compatible with the Suzaku observation in 2007 (Fig. 
[4}. The addition of a CS redwards of the Iron Ka line core 
is required by the flt (Ax^ ~ 20 with a signiflcance greater 
than 99.99%) and its fiux, being 24 ± 7% of the fiux of the 
narrow core of the Fe Ka, is consistent with expectations. 
Both Iron Ka and CS fluxes are in agree ment with the one s 
found in the previous XMM observation (|Matt et al.ll2004h . 
The equivalent widths of the high energy emission lines are 
in full agreement with the ones we flnd when we analyse the 
5 sets of data separately. We only find upper limits on fiuxes 
(< 0.5 X 10"^ ph cm-2 s-i) and EWs (< 50 eV) of the Fe 
XXVI emission lines at 6.966 keV. 

As a last cross-check, we let the 5 different absorbing 
column density parameters free to vary in our fit, to check 
whether a possible variation may have been occurred in the 
1.5 months monitoring. The best fit values of the 5 different 
Nh do not show any significant variation with respect to the 
best flt value for the whole data set (9.0 ± 0.5 x 10^^ cm"^) 
with a non signiflcant improvement of the flt (Ax^ ~ 5 with 
four more parameters and a significance lower than 8%). 
This result brings further evidence to the argument that 
absorption variability on short time scales (hours-days) can 
be ruled out in our analysis of NGC 4507. If the intervening 
absorbing material had varied on such short time scales we 
would not have measured a constant column density in a 1.5 
months monitoring. The measured values would have been 
completely scattered over the range of values observed in 
the past (4 - 9 X lO^^cm"^). 



^ The F-test is not a reliable test for the significance of emis- 
sion lines unless their normalizations are allowed to be negative 
dProtassov et al.ll20o3) . 



NGC 4507: 3-10 keV best fits 




Figure 2. 3-10 kcV best fits: each spectrum is divided by the 
effective area of the instrument. The impact of the variation in 
the absorbing column density on the three spectra can be clearly 
seen. 



3.2 Chandra ACIS spectral analysis 

The baseline model is the same we used to fit the XMM- 
Newton data. The overall fit is very good (x^ =327/329) and 
the addition of a further emission line at 1.81±0.02 keV is re- 
quired (Ax^ = 14 with a significance greater than 99.99%), 
with a fiux of 4^2 ^ 10~® phcm~^ s"'^, marginally consis- 
tent with the emission from Si Ka already found in XMM- 
Newton best fits. The soft X-ray spectrum is produced by 
two photoionised phases, while the contribution by a colli- 
sional gas is not required by the fit. Equivalent widths, fiuxes 
and energy centroids of the emission lines found in the 5 - 
7.5 keV energy range are fully consistent with the results 
reported above. Best fit values are shown in Table [T] and the 
column is labeled as Obs. 6. 

The refiected primary continuum has been fitted as de- 
scribed before (see Sect. 13. If : the best fit value of the ab- 
sorbing column density is 6.5 ± 0.7 x 10^"^ cm~'^, leading to 
a 2.5 X 10^'' cm~^ variation at a 3ct confldence level in a 
time scale ranging between 1.5 (time interval between the 
first and the fifth XMM-Newton observation) and 4 months 
(time interval between the last XMM-Newton observation 
and the one with Chandra). In Fig. [2] we show the infiuence 
of the change in Nh on the spectral shape between 3 and 10 
keV. Considering the partial degeneracy between the column 
density and the spectral slope, the significance of the vari- 
ation is even stronger, as illustrated by the F-Nj^ contour 
plots shown in Fig. 3. 



4 PHYSICAL DISCUSSION 

From the X-ray data analysis presented above a column den- 
sity variation in a time scale of months is evident. We are 
going to describe, in the following, the physical implications 
of this result. Changes in absorbing column density are due 
to two different physical processes: a variation of the ion- 
izing primary radiation, which causes the variation in the 
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Parameter 


Obs. 1 


Obs. 2 


Obs. 3 


Obs. 4 


Obs. 5 


Obs. 6 


Set 1 


K (PN-MOS) 






1 02+°-"^ 




0.93l|;°3 






A^H (10^3 cm-2) 




9.7t«:9 


7 6+^-" 
' -"-1.3 


9.4ti;l 


8.otg:« 


fi c;+0.7 
D.O_Q 7 


^-'-'-0.5 


r 


1 8+°-^ 
-^■°-o.i 


2 0+S? 
"—0.3 


1 4+°'! 


1 9 + n ^ 


1 6+S-? 


1 5+of 


^-°-0.2 


F (CSg.S kev) 


< 1.5 


< 1.7 


^•''-0.9 




o-8tr4 


2.0+'i';? 




Flux 
EW 


44+0.3 


f. 40+0.01 
6.4U_o 0^ 

n + 0.5 
■^■*'-0.5 


fi 41+0 02 
°-4^-0.02 

„+0.7 
' -0.7 

245^45 


fi 41+0-01 
O-4-'--0.02 
4 0+0.6 

285l« 


4n+0.02 

4 5+0.5 
335t« 


f. 40+0.02 
O-4U-0.02 
4 0+0-9 

220l|0 


fi 401+0-003 
0-40i_o 003 
4 2+0.3 


Fe K/3 E 


7.058 


7.058 


7.058 


7.058 


7.058 


7.058 


7.058 


Flux 


< 0.6 


< 0.5 


< 0.6 


< 0.5 


< 0.8 


< 1.1 


3+''-2 

U.J_o 2 


EW 


< 50 


< 35 


< 45 


< 40 


< 70 


< 70 


Fe XXV E 


6.700 


6.700 


6.700 


6.700 


6.700 


6.700 


6.700 


Flux 


6+''-3 


< 0.3 




< 0.7 


< 0.3 


< 0.4 


'-'-•^-0.1 


EW 


.^4-20 
4OI25 


< 20 


. -+25 
45I25 


< 45 


< 20 


< 20 




kT (keV) 


Q 43+0.12 

"■^■^-0.08 


53+°-^" 

— 0.25 


"■^^-0.15 


"•^■^-0.15 


^^•^^-0.12 




43+0.11 

"•^■^-0.07 


logfJi 


0.05 




1.75+°-?^ 

' ^ — 0.10 


1.59+°-" 

^"^ — 0.07 


1.67+?,-?^ 

^ ' —0.10 


, oc: + 0.20 
^^—0.12 


1.64+0-01 

^ — 0.04 


log Njji 




20 9+S i? 


20 9+S 

^^■^—0.2 


21 i+°1 

—0.3 


20 8+° t 
" "— 0.3 


21 9+°-^ 


21.0+Oi 
"-0.2 


log (72 


_Qgq+0.10 
u.oy^o 25 


n 1R+0.28 
-U.l»_o 23 


_Q 45+0.30 
'^•*^-0.25 


_0 47+0-64 
'-'-^'-0.82 


'^••^'J-0.50 




-0 15+0-21 
'-'-l''-0.25 


log Nh2 


19 9+°-i 


20 4+'''5 


< 19.8 


19 8+'''^ 


10 Q+0.3 
-^^•^-0.5 


21.5«/i 


10 Q+0-2 
^^•^-0.3 


Fo.5-2 koV 


^■^-0.1 


0.4t«- 






^■^-0.1 


0.3«-2 




F2-IO kcV 


7 7+0.3 
' ■ ' -0.3 


8.0l«:^ 




8.o«:^3 


7 c: + 0.8 


io.o«-i 




xVd.o.f. 


308/292 


347/268 


242/272 


267/267 


280/283 


313/327 


1576/1459 



Table 1. Best fit values. Energies are in keV, line fluxes in 10~^ pli cm~2 s~i, observed fluxes in 10~12 erg cm~2 g-l and EWs in eV. 
Photoionisation parameters log(7i, log{/2 and column densities logA^ni, logA'^H2 are the best fit values of the two photoionised phases 
needed to reproduce the soft emission; kT is the energy of the additional coUisional phase (see text for detail) . 

ionization state of the absorber or variations in the amount 
of absorbing gas along the line of sight. In the case of NGC 
4507 the first physical scenario can be clearly ruled out be- 
cause the difference in the 2-10 keV fluxes between the ob- 
servations is not significant enough to justify a variation in 
the primary ionizing radiation. 

The black hole mass of NGC 4507 is estimated by 
means of stellar velo city dispersion to be 4.5 x lO^M© 
iMarinucci et al.ll2012l V We assume the dimensions of the 
X-ray emitting source Ds to be 10 Rg- This is in agreement 
with continuum variability studies and disk-corona emission 
models, all suggesting a compact central X-ray source, 
confined within a few Rg from the central black hole. We 
also assume that the size of the obscuring cloud Dc^Ds- 
A schematic view of the geometrical structure is shown in 
Fig. [5l The transverse velocity vk for one obscuring cloud 
is then simply given by the linear dimension of the X-ray 
source, Ds, divided by the crossing time Tcr'- 



Dc 10 GMbfa _i 1 

VK ^ — ^ — ^ ^ 70 km s M7.65T4, , (2) 

where we introduced the adimensional parameters M7.65 = 
Mbh/IQ'-'^^MQ and T4 ~ 1 x lO^s = 4 months. 
If we then consider the absorbing material located at a dis- 
tance R from the central X-ray source, moving with Keple- 
rian velocity, we can calculate R with the simple formula: 

R = — — = 2 p2 - 40 pc M7.65-R10 T4 , (3) 

where Rio = Ds/W Rg- In the case of NGC 4507 the lower 
limit on the crossing time is 1.5 months, the time interval 
between the first and last observation of our XMM-Newton 
campaign, during which the column density is maximal and 
nearly constant. The upper limit to the uncovering time 
(~ Ds/vk) is 4 months, time interval between the fifth 
XMM-Newton observation (August 2011) and the Chandra 
one (December 2011). 
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Figure 3. F-Nh contour plots of the 5 combined XMM-Newton 
observations and tiie Chandra December 2010 observation. Solid 
black, red and green lines corresponds to Icr, 2cr, 3cr confidence 
levels, respectively. 



Using these limits in the relations above we get a lower 
limit oi R — 7 pc MtssRiq ^^'^ upper limit oi R — 
40 pc Mr.esRio ■ These distances imply that the obscur- 
ing material is located well outside the BLR. To be con- 
sistent with the BLR location the obscuring cloud should 
have a linear size Dc > 20Dg . If we assume a typical 
BLR density ric ~ 10® cm"^ (|Osterbrockl 1 19891 ) and a 
Nh = 2.5 X lO^^cm"^ (difference between XMM-Newton 
and Chandra spectra) we get a linear size for the obscuring 
cloud Dc < 3.5 Ds'- these two values are clearly inconsis- 
tent. The Suzaku observation is the only one that really 
probes the typical BLR timescales in NGC 4507, but with 
an average column density of ~ 10^* cm~^ it is not possible 



to disentangle any inner component with N^r ~ 10 



10^ 



We reduced and analyzed two lOks long Swift-XRT ob- 
servations on the 24**" and SO**" of December 2005 (obsid 
00035465001 and 00035465003 respectively), to check for 
column density variations on timescales of days. The base- 
line model we used to reproduce the data is a simple ab- 
sorbed power law, a Gaussian emission line and a further 
soft power law for spectral features below ~3 keV. Data 
quality does not allow for a clear investigation of the F-Nh 
parameter space, so we adopted a fixed value of F — 1.8. 
Only a marginal column density variation (at Icr confidence 
level) is found between the two Swift observations and the 
two measurements of Nh are consistent at the 90% con- 
fidence level (Fig. |4]). The lack of significant variability on 
short time scales does not imply that a BLR component does 
not exist at all, it suggests indeed a much more complex en- 
vironment of abso r bing structures, as already discussed in 
iNardini fc Risalitil l|201ll 'l for the dwarf Seyfert galaxy NGC 
4395. 

The obscuring clouds' velocities we measured (70-170 km/s 
A^T.es) are at least one order of magnitude smaller than the 
typical BLR clouds' velocities. The absorption variability is 
due to circumnuclear material whic h is located at di stances 
consistent with the putative torus (|Antonuccilll993h . Such 
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Figure 4. Column density light curves, (a) From 1990 until now 
the s ource has presented several A^h variations. Values are taken 
from iRisaliti et al.l ||2002| '). A clear changing-look on time scales 
of years has been already discussed in B12, while for the first 
time a column density variation has been observed on time scales 
of months, (b) Column density light curve from June 2010 to 
December 2010. The 5 XMM-Newton observations do not show 
any evidence of variation in Nh, while the Chandra value clearly 
differs from the XMM combined best fit value, which is plotted 
as a red solid line with errors as red dashed lines. 



material cannot be located at much larger scales (i.e. dust 
lanes) since the NLR is not significantly affect ed by redden- 
ing, a s inferred by the observed Ha/H/3 ratio (jKewlev et alJ 

iooj. 

Since the presence of a Compton-thick material around the 
nucleus is invariably accompanied by a neutral iron narrow 
Ka-emission line and a cold refiection emission it is inter- 
esting to point out the fact that the spatial scale of the 
obscuring material in NGC 4507 is consistent with the dis- 
tance of the reflector observed with Chandr a in the nearby 
Compton-thick Sy2 NGC 4945 jMarinucci et al.ll20l3 ). 
The existence of a more complex structure surrounding the 
central engine of AGN rather than the one predicted by the 
uniflcation model has been proposed and widely discussed in 
the past few years tMaiolino fc Riek e 1995: Elvis 2000: Mat3 
I2OO0I) a nd recently inlBianchi et al.1 (|2007() : iRisaliti fc ElvisI 
l|2010h : [Bianchi et al.1 (|2012[ ): lElvisri|2012h . 
Accordingly to these models, only absorbing matter on very 
different scales can be responsible for the wide phenomenol- 
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1|K, 



^ "Stands rdn. 
torys 



BLR 



terial is located well outside the BLR, and suggest a much 
more complex environment of absorbing structures; 

• the distances we inferred suggest that a single, univer- 
sal structure of the absorber is not enough to reproduce the 
X-ray absorption variability of this AGN. Different reflec- 
tors/absorbers are responsible for the observed X-ray fea- 
tures. 



In the next future, following the results presented in lRisalitil 
(2002*), a monitoring of all the sources that have shown 
changes on time scales of years but none on shorter (hours- 
days) can be performed. A broad band, time-resolved study 
of AGN is fundamental for a better understanding of the 
complex circumnuclear material and its interaction and re- 
sponse to the primary radiation. 



Figure 5. Schematic view of the circumnuclear absorbing struc- 
ture. 



ogy of column density variations, leading to an overall sce- 
nario where different absorbers/reflectors are responsible for 
the spectral changes in the Seyfert 2 galaxies observed so far. 
The column density variation we measured in NGC 4507 is 
a further piece that can be added to the puzzle. In our anal- 
ysis the lack of any variation on short time scales (hours, 
days) excludes an absorber located in the BLR while the 
change on time scales of months leads to an absorber much 
farther from the X-ray source, differently with respect to 
other variable objects (e . g. NGC 13 6 5, NGC 4151 UGC 
4203: iRisaliti et aLllioOsl . |2007| . l2009l : iPuccetti et alJlioOTl : 
iRisaliti et al.ll2010l ) on long time scales (months, years) and 
rapidly changing on short ones observed so far. Our analy- 
sis provides further evidence that a universal circumnuclear 
structure of absorbing matter is therefore not suited for tak- 
ing into account all the observed phenomenology on absorp- 
tion variability in AGN. While it is true that absorption 
must occur on different scales, not all the objects present 
evidence of absorption from all the possible scales. 



5 CONCLUSIONS 

We reported in this paper the analysis of 5 XMM-Newton 
observations spanning a period of 6 weeks and a Chandra 
observation performed 4 months afterwards of the obscured 
AGN in NGC 4507. This source had shown strong column 
density variations in time scales of years (from 1990 until 
2007) and none in shorter time scales during the 3 days 
of Suzaku monitoring. We therefore investigated time scales 
ranging from 1.5 up to 4 months, looking for absorbing struc- 
tures located farther from the innermost X-ray source. 
Our results can be summarized as follows: 

• a column density variation of |AA^h| = 2.5 x 10^^ 
cm~^ at a So" confidence level on a time interval between 
1.5 and 4 months has been measured. Such time scales 
lead to distances of the absorber ranging from R = {7 — 
40) M^ssRio P'^ with corresponding velocities of 70 — 170 
km/s Afy.es. These distances imply that the obscuring ma- 
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